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SUMMARY 



Methods are described for calculating the period and 
frequency of vibration of cantilever wings and similar 
structures in which the weight and moment of inertia vary 
along the span. - Both the beam and torsional frcq.ucn.cie8 
may "bo calculated by these methods. The procedure is 
illustrated by examples. 

It is shown that a surprisingly close approximation 
to the beam frequency may be obtained by a very brief 
calculation in which the curvature of the wing in.vibra-n 
tion is assumed to be constant. A somewhat longer compu- 
tation permits taking account of the true curvature of 
the beam by a series of successive approximations which 
arc shown to be strongly convergent. 

Analogous methods are applied to calculations of the 
torsional frequency. For the first approximation it is 
assumed that the angle of twist varies linearly along the 
senispan. The true variation of the twist is computed by 
successive approximations which arc strongly convergent, 
as in the case of beam vibrations. - - 





No tat ion 




v, 


strain energy. 




K, 


kinetic energy. 




E, 


modulus of elasticity in 


tension. 


c » 


modulus of elasticity in 


shear . 


g. 


acceleration of gravity. 
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period of vibration. 
' f requency of vibration, 
distance from root, 
dof 1 ec t ion . 
radius of curvature, 
velocity in beam vibration, 
moment of inertia (beam), 
woight per unit length. 

polar moment of inortla of weight per unit length, 
concentrated -woight. 
shearing force 
bending moment, 
■torsional moment, 
angl e of --twist . 

angular -velocity in torsional vibration, 
length of cantilever. 

area inclosed within cross-section of torque 
resistant structure. 

perimeter of cross section. 

thickness of shell. 



Beam Vibration 

At the instant of maximum amplitude of vibration, 
the beam has no motion and no kinetic energy. Tho strain 
energy in the bent beam is given by: 



T= (1/2) / EI(d 2 y/dx 3 ) s dx 



(1) 
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.The kinetic energy "at zero strain on.ergy is: 
T ' " K = ( l/ 2g) / wv 2 dx . __ . (2) 

"*"* In simple harmonic motion:' 

v = 2rry / T (3) 
Substituting this value of v in (2) gives i 

K = (2-rr2/gT 3 ) / wy 2 dx (4) 
Y = K (5) 

whence: 

4tt 3 / wy 3 dx , v 

T 3 = _ (6) 

gf El(d 3 y/dx 3 ) 3 dx 

E = 1/3? (7) 

v- The distributed forces on the vibrating wing are 

eq.ua! to the mass times the acceleration. The accelera- ' 

tion is proportional to the deflection, so that the forces 
are proportional to tho weight times the deflection. It is 

necessary only to determine the relative magnitudes of the 
dcfloctions along the "beam. Their absolute magnitudes do 
not influence the period. 

Tho true form of the deflection curve can bo computed 
quito accurately and readily by a series of successive 
approximations, beginning with a simple form of curve, such 
as one with a constant radius, i.e.: 

d^y/dxs = c ons taut 

In each successive approximation to tho defloction 
curve, it is assumed "that tho distributed forces are pro- 
portional to the distributed weights multiplied by tho de- 
flections computed in tho preceding approximation. "The 
series of curves obtained in this manner converges very 
rapidly toward a final form. The calculations may be 
readily carried out by tabulations without actually draw- 
ing any curves. On the other hand, a graphical solution 
with the assistance of a mechanical integraph may be 
quicker than tabulations. 
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The procedure is shown by tho- following example • 

Example . 1 . - Compute the period of vibration in the 
beam direction of a tapered cantilever wing in which the 
moments of inertia and the running woights at eleven 
oqually spaced stations are as in table I. (Station num- 
bers represent distances from the root.) 

E = 10 7 lb. /in. 3 



Method of Procedure 

For the first approximation, the doflection yj_ in 
the last column of tablo I is assumed to be proportional 
to the square of the distance from the root. 

The second approximation y 3 is calculated in table 

II by taking the running loads equal to wy lt The shear 
and bending aro then computod by successive integrations. 
Since relative and not actual valuos aro required, and 
the station spacings are constant, the integrations con- 
sist. merely of summations, as will bo apparent from in- 
spection of the tablo, without multiplying by the station 
intervals . 

In accordance with beam theory; 

d 3 y/dx s = M/EI 

The slope dy/dx is determined simply by summa- 
tion of M/l, omitting division by the constant E and 
multiplication by the constant intorval. The deflection 
y 2 is obtained by summation of dy/dx. 

Tho summations for dy/dx and y aro from tho root 
outward (i.e., from the bottom upward in tho table), whilo 
the summations for S and M are from tho tip inward, 
•or from the top downward in the table. 

The variations in the ratio Vs/Yx in the la9t column 
of tablo II show the differences botwoen the first and sec- 
ond approximations to tho form of tho defloction curve. 

The third approximation y 3 is determined in table 

III by tho sane procedure as in tablo II. The ratio 

is practically constant from tip to root, showing that the 
curve has settled down to its final form. 
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In table IV, d 3 y/dx 3 is the M/ 1 of ta"ble III* The 

calculation and summation of l(d s y/dx 3 ) 2 from d 2 y/dx 2 

aro obvious.. The calculation and summation of wy» from 
y of table III arc also obvious. 

Substituting in equation (6) the numerical values 
of the summations in table IV: 

T 2 = 4tt3 X 65517 X 16 X IP 7 
386 X To 7 X 1916OC00 

= 0.005595 sec. 2 
T = 0.0748 sec. 
N = 13.56 v p s 



Approximate Solution 

If it is assumed that d 2 y/dx 2 = 1, corresponding 
to a constant radius of curvaturo, the calculation of T 
becomes tho extremely simple and self-explanatory process 
shown in table V. 

l(d 2 y/dx 3 ) 2 = I 

dy/dx divided by, the s'tation interval equals l t 2, 3j 4, 
etc., from tho bottom up, and y divided by the square 
of tho interval is the summation of dy/dx from the bot- 
tom up.- 

T 3 - 4tt3 x 10067 X 160000 
386 X TO" 7 x 3965 

= 0.005556 sec . 3 ■ • . 

T = 0.0745 sec. 

N = 13.42 v p s 

"he error in T and N is only 0.4 percont. 

Exan-ple 2. - Since the very small error in the approxi- 
mate method might be purely chance "in this example, another 
calculation was made, using the same moments of inertia as 
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in the previous example, "but -with an entirely different 
distribution of weight* iPhe' calculat.i on of GJ by suc- 
cessive approximation is carried' out in tables .71 to X, 
inclusive, carrying the approximation to y one step 
further than "before. From the summations in table X, 

2 4tt 3 X 18218 X 16 x 10 7 

T = =r= 

386 X 10 7 X 1267700 

T 3 = 0., 02352 sec. 3 

T = 0.1534 see. 

1? = 6.52 v p s 

Computation of 3? on the assumption that d 2 y/dx 3 = 1 
is carried out in table XI. 

2 4tt 2 X 42576 X 160000 

T = — 

386 X 10. X 2965 

= 0.02350 sec. 3 

T = 0.1533 sec. 

N '= 6. 52 v p s 

The error by the approximate .solution is even les,s. 
than in the preceding example. 

From those two examples, it may bo concluded that 
tho approximate method is satisfactory for most cantilovor 
wings* It may bo noted from tables IV and X that in tho 
first example the true values of d 3 y/dx 3 vary about 80 
percent between the root and the tip, while in the second 
example the variation is only 50 percent, which accounts 
for the closer results in the second examplo than in tho 
first; but it is still rather, remarkable that tho approxi- 
mate method should bo so accurate for curvature that varios 
from 50 percent to 80 percent from the aseumod constant 
valuo • 

Example 5. - Compute the period Of vibration of a 
•weightless cantilever .of uniform I = 300 in. 4 , carrying 
a concentrated 1 o.ad. :o£- .1 , 000 pounds at 200 inches from tho 
root. :E = 10 7 :1b . /in--. 3 . 
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The exact solution is: 

T = */~~y/2.lZ. ' 
where y is the static deflection at the weight. 

tfL 3 1000 x 200 3 

y = _ — - = „ = 0*89 in. 

EI '900 X 10 7 

T = a/0. 89/3. 13 = 0.30B sec. 
Assuming d 2 y/dx 2 - l/R = const. 

At tip, y = L 3 /2E = 2000Q/E 

/wy 2 dx = 1000 X (20000/R) 3 
/EI(d 2 y/dx 2 ) 2 dx = To 7 x 300 x 200/R 2 

rp 3 - 4t t 2 * 1000 * (3Q000/R) 3 



386 X TO 7 X 300 X 200/E 2 

= 0.068 sec. 2 

T = 0.261 sec. 

In this case the error is large, as might he expected 
from the fact that the actual d 2 y/dx 3 varies as x instead 
of "being constant.. In other words, for a good approxima- 
tion to the period or frequency, the as sum e "37 curve - o"F "de- 
flection must not ho too far from roality, although as shown 
in the. examples it can he quite surprisingly far without 
approciable error in the result. 



Torsional Frequency 

The twist of a thin-walled closed section under tor- 
sional moments is given hy: 



46 V/"- < a > 



dx 4A 2 G 

J 

The strain energy of the twisted sections is given 

hy : 
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ay q ae 

dx 2 dx 



(9) 



Let . 

a 



k = (1/4A S ) / (l/t) ds (10) 
Combining (8), (9) , and (10) , 



V = - 



2 

The kinetic energy is given by: 

dK _ Pojf 27T 3 P8 g 
dx 2g gT 3 



k V dx J 



% - fir rpe3dI tl2> 



whence 



V = X 



T 2 = *r a /PB»ax ' (13) 



g& / (l/k) (d8/dx) <s dx 

The period may be calculated by a procedure analogous 
to that used for beam vibrations, beginning with the as- 
sumption that d8/dx = 1, and finding the true form of 
the twist by succossivo approximations. 

Example 4 . - Find the torsional period of vibration 
of the cantilever wing having the characteristics given 
in table XII. 

Let ..&-=» 4,0.00,000 lb. /it*- 2 

In the f-irst approximation, 8 is assumed propor- 
tional to the distance from the -root.. 

In the socond approximation, the applied torsional 
moments are taken oqual to P8- at each station. The 
total torquo Q 3 is obtained by summation of these mo- 
ments (table XIII). The twist is given by d8/dx = kq/0> 
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Summation of kQ, 2 from tho root outward gives 6 2 , omit- 
ting division "by G. 

It is seen from tho last column of tablo XIII that 
82/81 * s a reasonably constant ratio; "but the successive 
approximations are carried one stop farther in table XIV. 
The variation in 9 3 /9 2 is negligible. 

The summations of (l/k)(dQ/dx) n.nd P8 3 aro cal- 
culated in tablo XV, using the values of d8/dz and 6 
from table XIV. 

The numerical valuos from tablo XV aro substituted in 
equation (13), giving the torsional periodl 

a 4-rr 2 x 132464 X 4 X 10 5 

336 X 4 x 10 6 x 6559 x 10 3 

= 0.000206 sec. 3 

T = 0.01435 sec. 

II = 69.6 v p s 

Tho short-cut or approximate calculation "based on tho 
assumption that dG/dx =1 is carried out in tahle XVI. 
Only the figures in the last column and the summation of 
the second column requira any computation. All tho others 
may be written down directly. Substituting the numerical 
values in aquation (13), 

4-rr 2 X 82120 X 400 

I 3 = =15 . 

386 X 4 X 10 X 4157 

= 0.000202 sec. 2 
T = 0.0142 see. 
IT = 70.4 v p s 
Tho error is only 1.0 percent. 



Bureau of Aeronautics, iravy Department, 
Washington, D. C, February 1936. 
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